An adaptive biomolecular condensation response is conserved across environmentally divergent species.

Cells must sense and respond to sudden maladaptive environmental changes-stresses-to survive and thrive. Across eukaryotes, stresses such as heat shock trigger conserved responses: growth arrest, a specific transcriptional response, and biomolecular condensation of protein and mRNA into structures known as stress granules under severe stress. The composition, formation mechanism, adaptive significance, and even evolutionary conservation of these condensed structures remain enigmatic. Here we provide a remarkable view into stress-triggered condensation, its evolutionary conservation and tuning, and its integration into other well-studied aspects of the stress response. Using three morphologically near-identical budding yeast species adapted to different thermal environments and diverged by up to 100 million years, we show that proteome-scale biomolecular condensation is tuned to species-specific thermal niches, closely tracking corresponding growth and transcriptional responses. In each species, poly(A)-binding protein-a core marker of stress granules-condenses in isolation at species-specific temperatures, with conserved molecular features and conformational changes modulating condensation. From the ecological to the molecular scale, our results reveal previously unappreciated levels of evolutionary selection in the eukaryotic stress response, while establishing a rich, tractable system for further inquiry.

How hydrophobicity, side chains, and salt affect the dimensions of disordered proteins.

Despite the generally accepted role of the hydrophobic effect as the driving force for folding, many intrinsically disordered proteins (IDPs), including those with hydrophobic content typical of foldable proteins, behave nearly as self-avoiding random walks (SARWs) under physiological conditions. Here, we tested how temperature and ionic conditions influence the dimensions of the N-terminal domain of pertactin (PNt), an IDP with an amino acid composition typical of folded proteins. While PNt contracts somewhat with temperature, it nevertheless remains expanded over 10-58°C, with a Flory exponent, ν, >0.50. Both low and high ionic strength also produce contraction in PNt, but this contraction is mitigated by reducing charge segregation. With 46% glycine and low hydrophobicity, the reduced form of snow flea anti-freeze protein (red-sfAFP) is unaffected by temperature and ionic strength and persists as a near-SARW, ν ~ 0.54, arguing that the thermal contraction of PNt is due to stronger interactions between hydrophobic side chains. Additionally, red-sfAFP is a proxy for the polypeptide backbone, which has been thought to collapse in water. Increasing the glycine segregation in red-sfAFP had minimal effect on ν. Water remained a good solvent even with 21 consecutive glycine residues (ν > 0.5), and red-sfAFP variants lacked stable backbone hydrogen bonds according to hydrogen exchange. Similarly, changing glycine segregation has little impact on ν in other glycine-rich proteins. These findings underscore the generality that many disordered states can be expanded and unstructured, and that the hydrophobic effect alone is insufficient to drive significant chain collapse for typical protein sequences.

Transcriptome-wide mRNA condensation precedes stress granule formation and excludes stress-induced transcripts.

Stress-induced condensation of mRNA and proteins into stress granules is conserved across eukaryotes, yet the function, formation mechanisms, and relation to well-studied conserved transcriptional responses remain largely unresolved. Stress-induced exposure of ribosome-free mRNA following translational shutoff is thought to cause condensation by allowing new multivalent RNA-dependent interactions, with RNA length and associated interaction capacity driving increased condensation. Here we show that, in striking contrast, virtually all mRNA species condense in response to multiple unrelated stresses in budding yeast, length plays a minor role, and instead, stress-induced transcripts are preferentially excluded from condensates, enabling their selective translation. Using both endogenous genes and reporter constructs, we show that translation initiation blockade, rather than resulting ribosome-free RNA, causes condensation. These translation initiation-inhibited condensates (TIICs) are biochemically detectable even when stress granules, defined as microscopically visible foci, are absent or blocked. TIICs occur in unstressed yeast cells, and, during stress, grow before the appearance of visible stress granules. Stress-induced transcripts are excluded from TIICs primarily due to the timing of their expression, rather than their sequence features. Together, our results reveal a simple system by which cells redirect translational activity to newly synthesized transcripts during stress, with broad implications for cellular regulation in changing conditions.

HDX-MS finds that partial unfolding with sequential domain activation controls condensation of a cellular stress marker.

Eukaryotic cells form condensates to sense and adapt to their environment [S. F. Banani, H. O. Lee, A. A. Hyman, M. K. Rosen, Nat. Rev. Mol. Cell Biol. 18, 285-298 (2017), H. Yoo, C. Triandafillou, D. A. Drummond, J. Biol. Chem. 294, 7151-7159 (2019)]. Poly(A)-binding protein (Pab1), a canonical stress granule marker, condenses upon heat shock or starvation, promoting adaptation [J. A. Riback et al., Cell 168, 1028-1040.e19 (2017)]. The molecular basis of condensation has remained elusive due to a dearth of techniques to probe structure directly in condensates. We apply hydrogen-deuterium exchange/mass spectrometry to investigate the mechanism of Pab1's condensation. Pab1's four RNA recognition motifs (RRMs) undergo different levels of partial unfolding upon condensation, and the changes are similar for thermal and pH stresses. Although structural heterogeneity is observed, the ability of MS to describe populations allows us to identify which regions contribute to the condensate's interaction network. Our data yield a picture of Pab1's stress-triggered condensation, which we term sequential activation (Fig. 1A), wherein each RRM becomes activated at a temperature where it partially unfolds and associates with other likewise activated RRMs to form the condensate. Subsequent association is dictated more by the underlying free energy surface than specific interactions, an effect we refer to as thermodynamic specificity. Our study represents an advance for elucidating the interactions that drive condensation. Furthermore, our findings demonstrate how condensation can use thermodynamic specificity to perform an acute response to multiple stresses, a potentially general mechanism for stress-responsive proteins.

Folding of prestin's anion-binding site and the mechanism of outer hair cell electromotility.

Prestin responds to transmembrane voltage fluctuations by changing its cross-sectional area, a process underlying the electromotility of outer hair cells and cochlear amplification. Prestin belongs to the SLC26 family of anion transporters yet is the only member capable of displaying electromotility. Prestin's voltage-dependent conformational changes are driven by the putative displacement of residue R399 and a set of sparse charged residues within the transmembrane domain, following the binding of a Cl- anion at a conserved binding site formed by the amino termini of the TM3 and TM10 helices. However, a major conundrum arises as to how an anion that binds in proximity to a positive charge (R399), can promote the voltage sensitivity of prestin. Using hydrogen-deuterium exchange mass spectrometry, we find that prestin displays an unstable anion-binding site, where folding of the amino termini of TM3 and TM10 is coupled to Cl- binding. This event shortens the TM3-TM10 electrostatic gap, thereby connecting the two helices, resulting in reduced cross-sectional area. These folding events upon anion binding are absent in SLC26A9, a non-electromotile transporter closely related to prestin. Dynamics of prestin embedded in a lipid bilayer closely match that in detergent micelle, except for a destabilized lipid-facing helix TM6 that is critical to prestin's mechanical expansion. We observe helix fraying at prestin's anion-binding site but cooperative unfolding of multiple lipid-facing helices, features that may promote prestin's fast electromechanical rearrangements. These results highlight a novel role of the folding equilibrium of the anion-binding site, and help define prestin's unique voltage-sensing mechanism and electromotility.

An adaptive biomolecular condensation response is conserved across environmentally divergent species.

Cells must sense and respond to sudden maladaptive environmental changes-stresses-to survive and thrive. Across eukaryotes, stresses such as heat shock trigger conserved responses: growth arrest, a specific transcriptional response, and biomolecular condensation of protein and mRNA into structures known as stress granules under severe stress. The composition, formation mechanism, adaptive significance, and even evolutionary conservation of these condensed structures remain enigmatic. Here we provide an unprecedented view into stress-triggered condensation, its evolutionary conservation and tuning, and its integration into other well-studied aspects of the stress response. Using three morphologically near-identical budding yeast species adapted to different thermal environments and diverged by up to 100 million years, we show that proteome-scale biomolecular condensation is tuned to species-specific thermal niches, closely tracking corresponding growth and transcriptional responses. In each species, poly(A)-binding protein-a core marker of stress granules-condenses in isolation at species-specific temperatures, with conserved molecular features and conformational changes modulating condensation. From the ecological to the molecular scale, our results reveal previously unappreciated levels of evolutionary selection in the eukaryotic stress response, while establishing a rich, tractable system for further inquiry.

Folding of Prestin's Anion-Binding Site and the Mechanism of Outer Hair Cell Electromotility.

Prestin responds to transmembrane voltage fluctuations by changing its cross-sectional area, a process underlying the electromotility of outer hair cells and cochlear amplification. Prestin belongs to the SLC26 family of anion transporters yet is the only member capable of displaying electromotility. Prestin's voltage-dependent conformational changes are driven by the putative displacement of residue R399 and a set of sparse charged residues within the transmembrane domain, following the binding of a Cl - anion at a conserved binding site formed by amino termini of the TM3 and TM10 helices. However, a major conundrum arises as to how an anion that binds in proximity to a positive charge (R399), can promote the voltage sensitivity of prestin. Using hydrogen-deuterium exchange mass spectrometry, we find that prestin displays an unstable anion-binding site, where folding of the amino termini of TM3 and TM10 is coupled to Cl - binding. This event shortens the TM3-TM10 electrostatic gap, thereby connecting the two helices, resulting in reduced cross-sectional area. These folding events upon anion-binding are absent in SLC26A9, a non-electromotile transporter closely related to prestin. Dynamics of prestin embedded in a lipid bilayer closely match that in detergent micelle, except for a destabilized lipid-facing helix TM6 that is critical to prestin's mechanical expansion. We observe helix fraying at prestin's anion-binding site but cooperative unfolding of multiple lipid-facing helices, features that may promote prestin's fast electromechanical rearrangements. These results highlight a novel role of the folding equilibrium of the anion-binding site, and helps define prestin's unique voltage-sensing mechanism and electromotility.

Factors That Control the Force Needed to Unfold a Membrane Protein in Silico Depend on the Mode of Denaturation.

Single-molecule force spectroscopy methods, such as AFM and magnetic tweezers, have proved extremely beneficial in elucidating folding pathways for soluble and membrane proteins. To identify factors that determine the force rupture levels in force-induced membrane protein unfolding, we applied our near-atomic-level Upside molecular dynamics package to study the vertical and lateral pulling of bacteriorhodopsin (bR) and GlpG, respectively. With our algorithm, we were able to selectively alter the magnitudes of individual interaction terms and identify that, for vertical pulling, hydrogen bond strength had the strongest effect, whereas other non-bonded protein and membrane-protein interactions had only moderate influences, except for the extraction of the last helix where the membrane-protein interactions had a stronger influence. The up-down topology of the transmembrane helices caused helices to be pulled out as pairs. The rate-limiting rupture event often was the loss of H-bonds and the ejection of the first helix, which then propagated tension to the second helix, which rapidly exited the bilayer. The pulling of the charged linkers across the membrane had minimal influence, as did changing the bilayer thickness. For the lateral pulling of GlpG, the rate-limiting rupture corresponded to the separation of the helices within the membrane, with the H-bonds generally being broken only afterward. Beyond providing a detailed picture of the rupture events, our study emphasizes that the pulling mode greatly affects the factors that determine the forces needed to unfold a membrane protein.

LILAC: enhanced actin imaging with an optogenetic Lifeact.

Lifeact is a popular peptide-based label of actin filaments in live cells. We have designed an improved Lifeact variant, LILAC, that binds to actin in light using the LOV2 protein. Light control allows the user to modulate actin labeling, enabling image analysis that leverages modulation for an enhanced view of F-actin dynamics in cells. Furthermore, the tool reduces actin perturbations and cell sickness caused by Lifeact overexpression.

Development of in vivo HDX-MS with applications to a TonB-dependent transporter and other proteins.

Hydrogen-deuterium exchange mass spectrometry (HDX-MS) is a powerful tool that monitors protein dynamics in solution. However, the reversible nature of HDX labels has largely limited the application to in vitro systems. Here, we describe a protocol for measuring HDX-MS in living Escherichia coli cells applied to BtuB, a TonB-dependent transporter found in outer membranes (OMs). BtuB is a convenient and biologically interesting system for testing in vivo HDX-MS due to its controllable HDX behavior and large structural rearrangements that occur during the B12 transport cycle. Our previous HDX-MS study in native OMs provided evidence for B12 binding and breaking of a salt bridge termed the Ionic Lock, an event that leads to the unfolding of the amino terminus. Although purified OMs provide a more native-like environment than reconstituted systems, disruption of the cell envelope during lysis perturbs the linkage between BtuB and the TonB complex that drives B12 transport. The in vivo HDX response of BtuB's plug domain (BtuBp) to B12 binding corroborates our previous in vitro findings that B12 alone is sufficient to break the Ionic Lock. In addition, we still find no evidence of B12 binding-induced unfolding in other regions of BtuBp that could enable B12 passage. Our protocol was successful in reporting on the HDX of several endogenous E. coli proteins measured in the same measurement. Our success in performing HDX in live cells opens the possibility for future HDX-MS studies in a native cellular environment. IMPORTANCE: We present a protocol for performing in vivo HDX-MS, focusing on BtuB, a protein whose native membrane environment is believed to be mechanistically important for B12 transport. The in vivo HDX-MS data corroborate the conclusions from our previous in vitro HDX-MS study of the allostery initiated by B12 binding. Our success with BtuB and other proteins opens the possibility for performing additional HDX-MS studies in a native cellular environment.

HDX-MS performed on BtuB in E. coli outer membranes delineates the luminal domain's allostery and unfolding upon B12 and TonB binding.

To import large metabolites across the outer membrane of gram-negative bacteria, TonB-dependent transporters (TBDTs) undergo significant conformational change. After substrate binding in BtuB, the Escherichia coli vitamin B12 TBDT, TonB binds and couples BtuB to the inner-membrane proton motive force that powers transport [N. Noinaj, M. Guillier, T. J. Barnard, S. K. Buchanan, Annu. Rev. Microbiol. 64, 43­60 (2010)]. However, the role of TonB in rearranging the plug domain of BtuB to form a putative pore remains enigmatic. Some studies focus on force-mediated unfolding [S. J. Hickman, R. E. M. Cooper, L. Bellucci, E. Paci, D. J. Brockwell, Nat. Commun. 8, 14804 (2017)], while others propose force-independent pore formation by TonB binding [T. D. Nilaweera, D. A. Nyenhuis, D. S. Cafiso, eLife 10, e68548 (2021)], leading to breakage of a salt bridge termed the "Ionic Lock." Our hydrogen­deuterium exchange/mass spectrometry (HDX-MS) measurements in E. coli outer membranes find that the region surrounding the Ionic Lock, far from the B12 site, is fully destabilized upon substrate binding. A comparison of the exchange between the B12-bound and the B12+TonB­bound complexes indicates that B12 binding is sufficient to unfold the Ionic Lock region, with the subsequent binding of a TonB fragment having much weaker effects. TonB binding accelerates exchange in the third substrate-binding loop, but pore formation does not obviously occur in this or any region. This study provides a detailed structural and energetic description of the early stages of B12 passage that provides support both for and against current models of the transport process.

Challenges and Advantages of Accounting for Backbone Flexibility in Prediction of Protein-Protein Complexes.

Predicting protein binding is a core problem of computational biophysics. That this objective can be partly achieved with some amount of success using docking algorithms based on rigid protein models is remarkable, although going further requires allowing for protein flexibility. However, accurately capturing the conformational changes upon binding remains an enduring challenge for docking algorithms. Here, we adapt our Upside folding model, where side chains are represented as multi-position beads, to explore how flexibility may impact predictions of protein-protein complexes. Specifically, the Upside model is used to investigate where backbone flexibility helps, which types of interactions are important, and what is the impact of coarse graining. These efforts also shed light on the relative challenges posed by folding and docking. After training the Upside energy function for docking, the model is competitive with the established all-atom methods. However, allowing for backbone flexibility during docking is generally detrimental, as the presence of comparatively minor (3-5 Å) deviations relative to the docked structure has a large negative effect on performance. While this issue appears to be inherent to current forcefield-guided flexible docking methods, systems involving the co-folding of flexible loops such as antibody-antigen complexes represent an interesting exception. In this case, binding is improved when backbone flexibility is allowed using the Upside model.

Prediction and Validation of a Protein's Free Energy Surface Using Hydrogen Exchange and (Importantly) Its Denaturant Dependence.

The denaturant dependence of hydrogen-deuterium exchange (HDX) is a powerful measurement to identify the breaking of individual H-bonds and map the free energy surface (FES) of a protein including the very rare states. Molecular dynamics (MD) can identify each partial unfolding event with atomic-level resolution. Hence, their combination provides a great opportunity to test the accuracy of simulations and to verify the interpretation of HDX data. For this comparison, we use Upside, our new and extremely fast MD package that is capable of folding proteins with an accuracy comparable to that of all-atom methods. The FESs of two naturally occurring and two designed proteins are so generated and compared to our NMR/HDX data. We find that Upside's accuracy is considerably improved upon modifying the energy function using a new machine-learning procedure that trains for proper protein behavior including realistic denatured states in addition to stable native states. The resulting increase in cooperativity is critical for replicating the HDX data and protein stability, indicating that we have properly encoded the underlying physiochemical interactions into an MD package. We did observe some mismatch, however, underscoring the ongoing challenges faced by simulations in calculating accurate FESs. Nevertheless, our ensembles can identify the properties of the fluctuations that lead to HDX, whether they be small-, medium-, or large-scale openings, and can speak to the breadth of the native ensemble that has been a matter of debate.

Lipid bilayer induces contraction of the denatured state ensemble of a helical-bundle membrane protein.

Defining the denatured state ensemble (DSE) and disordered proteins is essential to understanding folding, chaperone action, degradation, and translocation. As compared with water-soluble proteins, the DSE of membrane proteins is much less characterized. Here, we measure the DSE of the helical membrane protein GlpG of Escherichia coli (E. coli) in native-like lipid bilayers. The DSE was obtained using our steric trapping method, which couples denaturation of doubly biotinylated GlpG to binding of two streptavidin molecules. The helices and loops are probed using limited proteolysis and mass spectrometry, while the dimensions are determined using our paramagnetic biotin derivative and double electron-electron resonance spectroscopy. These data, along with our Upside simulations, identify the DSE as being highly dynamic, involving the topology changes and unfolding of some of the transmembrane (TM) helices. The DSE is expanded relative to the native state but only to 15 to 75% of the fully expanded condition. The degree of expansion depends on the local protein packing and the lipid composition. E. coli's lipid bilayer promotes the association of TM helices in the DSE and, probably in general, facilitates interhelical interactions. This tendency may be the outcome of a general lipophobic effect of proteins within the cell membranes.

Engineered Metal-Binding Sites to Probe Protein Folding Transition States: Psi Analysis.

The formation of the transition state ensemble (TSE) represents the rate-limiting step in protein folding. The TSE is the least populated state on the pathway, and its characterization remains a challenge. Properties of the TSE can be inferred from the effects on folding and unfolding rates for various perturbations. A difficulty remains on how to translate these kinetic effects to structural properties of the TSE. Several factors can obscure the translation of point mutations in the frequently used method, "mutational Phi analysis." We take a complementary approach in "Psi analysis," employing rationally inserted metal binding sites designed to probe pairwise contacts in the TSE. These contacts can be confidently identified and used to construct structural models of the TSE. The method has been applied to multiple proteins and consistently produces a considerably more structured and native-like TSE than Phi analysis. This difference has significant implications to our understanding of protein folding mechanisms. Here we describe the application of the method and discuss how it can be used to study other conformational transitions such as binding.

Folding and misfolding of potassium channel monomers during assembly and tetramerization.

The dynamics and folding of potassium channel pore domain monomers are connected to the kinetics of tetramer assembly. In all-atom molecular dynamics simulations of Kv1.2 and KcsA channels, monomers adopt multiple nonnative conformations while the three helices remain folded. Consistent with this picture, NMR studies also find the monomers to be dynamic and structurally heterogeneous. However, a KcsA construct with a disulfide bridge engineered between the two transmembrane helices has an NMR spectrum with well-dispersed peaks, suggesting that the monomer can be locked into a native-like conformation that is similar to that observed in the folded tetramer. During tetramerization, fluoresence resonance energy transfer (FRET) data indicate that monomers rapidly oligomerize upon insertion into liposomes, likely forming a protein-dense region. Folding within this region occurs along separate fast and slow routes, with τfold ∼40 and 1,500 s, respectively. In contrast, constructs bearing the disulfide bond mainly fold via the faster pathway, suggesting that maintaining the transmembrane helices in their native orientation reduces misfolding. Interestingly, folding is concentration independent despite the tetrameric nature of the channel, indicating that the rate-limiting step is unimolecular and occurs after monomer association in the protein-dense region. We propose that the rapid formation of protein-dense regions may help with the assembly of multimeric membrane proteins by bringing together the nascent components prior to assembly. Finally, despite its name, the addition of KcsA's C-terminal "tetramerization" domain does not hasten the kinetics of tetramerization.

Properties of protein unfolded states suggest broad selection for expanded conformational ensembles.

Much attention is being paid to conformational biases in the ensembles of intrinsically disordered proteins. However, it is currently unknown whether or how conformational biases within the disordered ensembles of foldable proteins affect function in vivo. Recently, we demonstrated that water can be a good solvent for unfolded polypeptide chains, even those with a hydrophobic and charged sequence composition typical of folded proteins. These results run counter to the generally accepted model that protein folding begins with hydrophobicity-driven chain collapse. Here we investigate what other features, beyond amino acid composition, govern chain collapse. We found that local clustering of hydrophobic and/or charged residues leads to significant collapse of the unfolded ensemble of pertactin, a secreted autotransporter virulence protein from Bordetella pertussis, as measured by small angle X-ray scattering (SAXS). Sequence patterns that lead to collapse also correlate with increased intermolecular polypeptide chain association and aggregation. Crucially, sequence patterns that support an expanded conformational ensemble enhance pertactin secretion to the bacterial cell surface. Similar sequence pattern features are enriched across the large and diverse family of autotransporter virulence proteins, suggesting sequence patterns that favor an expanded conformational ensemble are under selection for efficient autotransporter protein secretion, a necessary prerequisite for virulence. More broadly, we found that sequence patterns that lead to more expanded conformational ensembles are enriched across water-soluble proteins in general, suggesting protein sequences are under selection to regulate collapse and minimize protein aggregation, in addition to their roles in stabilizing folded protein structures.

Water as a Good Solvent for Unfolded Proteins: Folding and Collapse are Fundamentally Different.

The argument that the hydrophobic effect is the primary effect driving the folding of globular proteins is nearly universally accepted (including by the authors). But does this view also imply that water is a "poor" solvent for the unfolded states of these same proteins? Here we argue that the answer is "no," that is, folding to a well-packed, extensively hydrogen-bonded native structure differs fundamentally from the nonspecific chain collapse that defines a poor solvent. Thus, the observation that a protein folds in water does not necessitate that water is a poor solvent for its unfolded state. Indeed, chain-solvent interactions that are marginally more favorable than nonspecific intrachain interactions are beneficial to protein function because they destabilize deleterious misfolded conformations and inter-chain interactions.

Structural basis for adhesion G protein-coupled receptor Gpr126 function.

Many drugs target the extracellular regions (ECRs) of cell-surface receptors. The large and alternatively-spliced ECRs of adhesion G protein-coupled receptors (aGPCRs) have key functions in diverse biological processes including neurodevelopment, embryogenesis, and tumorigenesis. However, their structures and mechanisms of action remain unclear, hampering drug development. The aGPCR Gpr126/Adgrg6 regulates Schwann cell myelination, ear canal formation, and heart development; and GPR126 mutations cause myelination defects in human. Here, we determine the structure of the complete zebrafish Gpr126 ECR and reveal five domains including a previously unknown domain. Strikingly, the Gpr126 ECR adopts a closed conformation that is stabilized by an alternatively spliced linker and a conserved calcium-binding site. Alternative splicing regulates ECR conformation and receptor signaling, while mutagenesis of the calcium-binding site abolishes Gpr126 function in vivo. These results demonstrate that Gpr126 ECR utilizes a multi-faceted dynamic approach to regulate receptor function and provide relevant insights for ECR-targeted drug design.